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ABSTRACT 

We present a study of the central 2.3 kpc of NGC 5135, a nearby Luminous Infrared Galaxy (LIRG) 
with an AGN and circumnuclear starburst. Our main results are based on intermediate spectral 
resolution 3000-4000) near infrared data taken with the SINFONI integral field spectrograph at 
the ESO VLT. The ionization of the different phases of the interstellar gas and the complex structures 
of the star formation have been mapped. Individual regions of interest have been identified and studied 
in detail. 

For the first time in this galaxy, we have detected the presence of a high excitation ionization cone 
centered on the AGN by using the [SiVI]A1.96 /xm line. So far, this structure is the largest reported in 
the literature for this coronal line, extending (in projection) as far as ~ 600 pc from the galaxy nucleus. 
In a complex spatial distribution, a variety of mechanisms are driving the gas ionization, including SNe 
remnant shocks, young stars and AGN photoionization. The excitation of the molecular gas, however, 
is mainly produced by X-rays and SNe remnant shocks. UV-mechanisms like fluorescence represent a 
marginal overall contribution to this process, contrary to the expectations we might have for a galaxy 
with a recent and strong star formation. Our SNe rate estimations from [FeII]A1.64 fj,ui are in excellent 
agreement with 6 cm radio emission predictions. Typical SNe rates between 0.01-0.04 yr~^ were found 
for individual « 200 pc-scale regions, with an overall SNe rate of 0.4-0.5 yr~^. Even though NGC 5135 
has suffered a recent starburst (6-7 Myr ago), the data strongly suggest the presence of a second, 
older stellar population dominated by red giant/supergiant stars. However, simple stellar population 
models cannot sharply discriminate between the different populations. 

Subject headings: galaxies: Seyfert - galaxies: nuclei - galaxies: structure - galaxies: ISM - infrared: 
galaxies - infrared: ISM - (ISM:) supernova remnants - galaxies: general 



1. INTRODUCTION 



Since their discovery (Kleinmann & Low 1970, Rieke & 
Low 1972), the importance of low-z Luminous {10^^ L(. 



< 



LiR < IO^^Lq, LIRG) and Ultraluminous {Lir > 
lO^^i©, ULIRG) Infrared Galaxies has been widely rec- 
ognized. (U)LIRGs are gas rich galaxies (Solomon et al. 
1997), where an active nucleus and powerful starbursts 
coexist and contribute to the energy output (e.g. Gen- 
zel et al. 1998, Spoon et al. 2007). While LIRGs ap- 
pear to be mostly spirals (Arribas et al. 2004, Alonso- 
Herrero et al. 2006a), ULIRGs are strongly interacting 
systems and mergers (e.g. Bushouse et al. 2002) evolv- 
ing into intermediate-mass ellipticals (e.g. Genzel et al. 
2001, Dasyra et al. 2006). Local (U)LIRGs have been 
proposed as possible counterparts of the submillimeter 
population observed at higher redshifts (Blain et al. 2002 
for a review, see also Chapman et al. 2003, Frayer et 
al. 2003, Egami et al. 2004). Also, cosmological sur- 
veys with Spitzer have shown that the majority of IR 
selected galaxies at z < 1 are in the LIRG class, while 
LIRGs and ULIRGs make a significant contribution to 
the IR galaxy population and to the star formation at 
1 < z < 2 and z > 2, respectively (Egami et al. 2004; 
Le Floch et al. 2004, 2005; Perez-Gonzalez et al. 2005; 
Caputi et al. 2007). 

Detailed investigations of the physical properties, stel- 
lar populations, AGN-starburst connection and gas flows 
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on these complex systems can only be obtained through 
integral field spectroscopy (IFS). Initial studies of small 
samples of (U)LIRGs based on 4-meter class telescope 
optical IFS have already been obtained (e.g. Colina et 
al. 2005, Monreal-Ibero et al. 2006, Alonso-Herrero et al. 
2009). To extend these studies to larger samples, and 
also to the near-IR, we have started a survey of low-z 
(U)LIRGs using state-of-the-art IFS like VLT/VIMOS 
(Arribas et al. 2008) and VLT/SINFONI (visual and 
near-IR, respectively). This survey (SIRIUS: Survey of 
luminous IR galaxies with Jntegral field Units) will allow 
us to characterize the kpc-scale ionization and kinematics 
of a representative sample of low-z (U)LIRGs covering a 
wide luminosity range, several morphologies from spirals 
to interacting and advanced mergers, as well as different 
classes of activity. This will also form a local reference 
for future IFS studies of high-z IR galaxies with instru- 
ments such as the Near-IR Spectrograph {NIRSpec) and 
Mid-IR Instrument [MIRI) on board of the James Webb 
Space Telescope (Gardner et al. 2006). 

As part of our survey, we present our first results 
with SINFONI showing the power of near-IR integral 
field spectroscopy by studying the local LIRG NGC 5135. 
This is an SBab galaxy at z = 0.01396 (from NED^, at 
w 58.7Mpc assuming Hq — 70kms^^ Mpc~^) which be- 
longs to a group of seven galaxies (KoUatschny & Fricke 
1989). Its dual nature as a starburst hosting an AGN 
and its almost face-on sky orientation make NGC 5135 
an ideal prototype-object for detailed studies of these 

^ http:/ /ncdwww. ipac.caltech.edu/ 
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hybrid systems. 

The first attempts to study NGC5135 was made 25 
years ago by Huchra(1983), Huchra et al. (1983) and 
PhiUips et al. (1983), who classified this object as a 
Seyfert2. A year later, Thuan (1984) found that this 
galaxy is also a bright UV source, with mixed Seyfert 
and starburst spectral characteristics (see also Kinney et 
al. 1993). Ulvestad & Wilson (1989) observed NGC5135 
at 6 and 20 cm using the VLA, as part of a larger sam- 
ple of nearby Seyfert galaxies. An important diffuse 
radio component was measured in this morphologically 
asymmetric object, which is in turn aligned with the ex- 
tended Ha emission (Haniff et al. 1988, Garcia Barreto et 
al. 1996). Because of the lack of an accurate optical nu- 
cleus, Ulvestad & Wilson chose the brightest radio source 
as the reference for the central region. 

Almost 10 years later, Gonzalez Delgado et al. (1998) 
published a study based on Hubble Space Telescope 
{HST from now on) and ground-based UV/optical imag- 
ing and spectroscopy to study possible starburst- AGN 
connections. They found independent pieces of evidence 
pointing towards the presence of a dusty, compact and 
powerful starburst in the central region (inner 2") of 
NGC5135, which dominates the UV spectrum. From 
spectral modeling and colors, these authors estimated 
an age of between 3 and 6 Myr for the last starburst 
episode. They also concluded that the warm interstellar 
gas is outflowing, most likely because of the starbursts. 
In the optical, these authors detected a nuclear bar for 
this galaxy, which had been previously found in the near 
infrared (hereafter NIR) by Mulchaey et al. (1997). 

Using Chandra X-ray observations, Levenson et al. 
(2004) revisited NGC 5135. Two strong central sources 
were detected and the (AGN) nucleus was correctly iden- 
tified for the first time, showing that its important in- 
ternal extinction had made it difhcult to recognize in 
previous works. They highlighted the importance of the 
stellar processes even in X-rays, showing that almost all 
the soft emission (0.4-1 keV) has a stellar origin, even 
nearby the AGN. Only the hard emission (4-10 keV) is 
dominated by the active nucleus. 

Alonso-Herrero et al. (2006a, b) and Dfaz-Santos et al. 
(2008) studied this object as part of a sample of ~ 30 
local LIRGs. With high resolution mid and NIR im- 
ages (Gemini, HST) in different bands, they identified 
the individual Hii knots which are hidden at shorter 
wavelengths. In NGC 5135, compact Paa emission (^ 1- 
2 kpc) from star formation was detected as well as a 
strong, nuclear warm dust component in the A'^band. 

The NGC 5135 SINFONI data allowed us to address 
some of the issues described above, as well as to show 
the scientific potential of our SINFONI data set. On 
the other hand, issues such as the general role of 
(U)LIRGs in galaxy formation and evolution will be 
better addressed with our entire sample in future papers 
of the series. We leave for the following paper of the 
series (hereafter. Paper II) a detailed kinematical study 
and mass estimation for this galaxy. 

The paper is organized as follows: in Sec. [2] the obser- 
vations ()2.1|) , data reduction ()2.2|1 and analysis (12. 3|) are 
described. Also, several maps are presented (|2.4|) and 
specific regions of study are identified (|2.5p . In Sec. [3] the 
main results and discussion arc shown, including a study 



of the central gas and stellar structure (|3.ip . internal ex- 
tinction structure (|3.2p . gas ionization mechanisms and 
young stellar population age (|3.3p . further characteriza- 
tion of the stellar population ages by using stellar absorp- 
tion features ()3.4|) . the II2 excitation mechanisms p.Sp 
and estimations of the supernova rate (|3.6p . Finally, in 
Sec. 21 we present the conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observing with SINFONI 

The observations were carried out in service mode 
with SINFONI (Eisenhauer et al. 2003, Bonnet et al. 
2004), the NIR integral field spectrometer mounted at 
VLT-Yepun (Cerro Paranal, Chile), with no AO-assisted 
guiding. The galaxy was observed separately in the H 
(1.45 - 1.85^m ) and K (1.95 - 2.45^m) bands during 
April and May 2006. Observing the H and if-bands 
independently provides a higher spectral resolution (~ 
3000 and ~ 4000, respectively) compared to the com- 
bined H+K ha,nd mode (- 1500). The full-width-at-half- 
maximum (FWHM from now on) as measured from sky 
lines is 6 A for the _ff-band and about ~ 5 A for K. The 
dispersions are 2 and 2.5 A pixel" ^ for 77 and K, respec- 
tively. The 64 X 64 spatial pixel (spaxels) field-of-view 
(FoV from now on) was centered close to the region where 
the AGN is located (Levenson et al. 2004). With a scale 
of 0.250"pixel~^, SINFONI provides a two-dimensional 
field of 8" X 8". 

For the i?-band galaxy data, the total on-target inte- 
gration time was 40 minutes (divided into 3 observation 
blocks, hereafter OBs), with an average seeing (FWHM 
of the point-spread-function, hereafter PSF) and airmass 
of 0.54" and 1.24, respectively. For the i^T-band data, 
the corresponding values were 25 minutes (2 OBs), 0.62" 
and 1.04 airmass. An entire iiT-band OB was discarded 
because of the high seeing (> 2"). Given the strong vari- 
ations of the IR sky with time, the observing blocks were 
split into shorter on-source integrations of 150 seconds 
each. Between them, we obtained identical exposures 
of the sky background (ABBA pattern) by nodding the 
telescope to a point shifted Aa = 42" and Ai5 = 87" 
from the galaxy nucleus. In this way, each "object in- 
tegration" has a corresponding "sky integration" which 
was subtracted during the pipeline processing to remove 
(to first order) the sky signal from the final data cube. 
Also, the object and sky frames were jittered within a box 
width of 1.5" to get rid of detector effects (bad pixels). 

2.2. Data Reduction 

The basic data reduction was performed with the SIN- 
FONI Pipeline (version 1.8.2). Dark subtraction, flat- 
field normalisation, pixel linearity, geometrical distortion 
and wavelength calibration corrections were applied for 
each object-sky pair by using the products of SINFONI 
standard calibration plan. The calibration products (e. g. 
darks and flats) are generated on a daily basis, with the 
exception of the linearity and geometrical distortion data 
which are produced monthly. The wavelength calibration 
was performed with Xe-Ar and Ne lamps for the H and 
A'-bands, respectively. The arc-frames include about 15 
strong emission lines, well spread along the correspond- 
ing wavelength ranges. We double-checked the wave- 
length calibrations on the individual object-sky frames 
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by tracing the position of about 10-15 strong sky-Unes 
along the detector. For the fl-band data, we found a 
systematic offset in wavelength of +1.6 ± O.lA, indepen- 
dent of the data frame used for the measurements. A 
similar situation affects the JT-band data, where an off- 
set of -t-2.3 ± 0.3A was found. Although these shifts are 
smaller than the spectral resolution (FWHM = 5-6A) 
they correspond to almost 1 pixel offset. Therefore, we 
decided to correct for them given their amplitude and 
systematic presence in the data. As a final check and 
after applying the offset corrections, we compared for 
several spaxels the measured and theoretical arc-line po- 
sitions, repeating the measurements in different frames. 
From the distribution of their differences, we estimated 
the (3 a) uncertainty in the wavelength calibration to be 
< 1.3 and < 1 A for the H- and K-hand spectral ranges, 
respectively. 

The background sky was subtracted for each individ- 
ual 150 s object exposure. The SINFONI Pipeline does 
not only subtract the object/sky pairs but it also applies 
other correction algorithms in order to minimise the sky 
line residuals. Despite the fact that the sky line relative 
positions of each object/sky pair were between spectral 
resolution limits, the sky frames were slightly shifted in 
wavelength with respect to their object pairs to reduce 
the P-Cygni-type residuals. In a second step, the pipeline 
includes an implementation of the sky subtraction pack- 
age described by Davies et al. (2007). This code allows 
to separately scale different transitions of the telluric OH 
lines before sky subtraction. After the correction, the 
sky-subtracted object frames present very small sky-line 
residuals whose importance (in the worst cases) is purely 
cosmetic. 

At this point, the signal (counts) from the object 
frames within an OB (between 5 and 6 on-target frames 
per OB) were combined to generate the integrated data- 
cubes for each band. Given the jittering applied dur- 
ing observations, the pipeline re-aligned the individual 
object frame centers before combination. As a test for 
these alignments, we generate a individual data-cubes 
from each object frame; the position of the three more 
prominent peaks in the continuum (for H and if-band) 
were traced by fitting two-dimensional elliptical Gaus- 
sians and the offsets were re-calculated. The data-cubes 
generated with the new offsets were in excellent agree- 
ment with the pipeline products. In fact, the majority 
of the offsets used by the pipeline were recovered by our 
method within a 10% uncertainty. Consequently, we de- 
cided to stick to the pipeline data-cube products, while 
our cube alignment method will be used in following steps 
of the reduction. 

As part of the calibration plan, spectrophotometric 
standard stars were observed along with our target ob- 
ject in order to perform the flux calibration. For each 
of the three iJ-band OBs an early-type star was ob- 
served (Hip 075902, Hip 068496 and Hip 065630 of spec- 
tral types B9V, B7II and B3IV, respectively), while 
only one star was used for the two consecutive JT-band 
OBs (Hip 078004, B3/B4V). The stellar spectra were 
extracted within an aperture of 6 a of the best two- 
dimensional Gaussian fit of the light profile (between 5 
and 6"). The measured spectra were divided by a theo- 
retical spectrum for the corresponding stellar type, gen- 



erating a multiplicative "sensitivity function" to be ap- 
phed to the object data. We use the models of Pickles 
(1998) as our theoretical templates after scaling them 
to the corresponding magnitudes listed in the Hipparcos 
catalogue (Ferryman et al. 1997). 

For the //-band stars, we found some differences be- 
tween the three individual sensitivity functions. The 
corrections in the central range 1.52-1.75 /im are in very 
good agreement between the three functions, showing 
relative deviations within 5%. In the extremes of the 
wavelength range (< 1.50 and > 1.75 /im), however, rel- 
ative differences could be as large as 15%. This sug- 
gests that the main source of uncertainty is not coming 
from Pickles' models (our "real" spectra) neither from 
peculiarities of the individual stars but from variations 
in the atmospheric transmission during observation (each 
star was observed just after one object OB, with airmass 
of 1.1, 1.3 and 1.5, respectively). We decided to apply 
the individual sensitivity functions to their correspond- 
ing OBs. Each of these functions has implicit the atmo- 
spheric transmission at a time and airmass close to that 
of the object observation. The results for the flux cali- 
brated OB's data were very good, with no strong resid- 
uals from transmission variations. In any case, because 
the only _H-band line we will use is the [FeII]A1.64 /im 
(essentially in the center of the band), the effects of the 
correction in the extremes of the wavelength range are 
purely cosmetic. 

One problem we found with Pickles' models is that for 
many of our spectral types the models do not include 
particular stellar features in the NIR spectrum, but they 
just present a smooth continuum resembling that of a 
black body. If we do not take care of this, the spectral 
features of the star will show up in the galaxy spectrum. 
Fortunately, lines from the Brackett series in absorption 
are essentially the only spectral features we found in our 
stellar spectra. So they were modeled with Gaussians 
and the stellar spectra were normalised before flux cali- 
bration. The residuals from this process were negligible. 

Once the OB's spectra were reduced and flux cali- 
brated, we combined the individual OBs in "final" data 
cubes for H- and i^T-bands. We estimated the offsets be- 
tween OBs by finding the three main peaks in the contin- 
uum, as described above for the individual frame combi- 
nation within an OB. As a final check, we compared our 
integrated fluxes in H and K with 2MASS data (Skrut- 
skie et al. 2006). After adapting the spectral range to 
the Ks-hand used in 2MASS, selecting a galaxy's central 
region comparable to our FoV and taking into account 
zero-point considerations, we found a good agreement 
in fluxes within 10% uncertainty for both bands. This 
is consistent with our maximum error estimation (15%) 
found between individual OB's flux calibrations in H. In 
Fig.[T]and[2]we present K- and iJ-band integrated spec- 
tra of six regions in NGC 5135 (see Fig. [3] and Sees. 12. 51 
and l3.1l for further details about the regions). 

2.3. Emission Line Fitting 

The (redshifted) line peak, flux, equivalent-width 
(EW) and FWHM were calculated from many of the 
emission lines presented in H- and if-bands. Other pa- 
rameters like galactocentric velocities and velocity dis- 
persions will be presented and analysed in Paper II. 

A single Gaussian component provides a good fit for 
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Fig. 1. — Observed-frame, normalized K-hand spectra of six selected regions in NGC5135. The spectra were integrated within 
circular apertures of 0.62" diameter, corresponding to ~ 180 pc. Different emission lines and absorption bands are labeled. Blue 
labeled features are those used in this study. 
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Fig. 2. — Observed- frame, normalized H-hand spectra of six selected regions in NGC 5135. The spectra were integrated within 
circular apertures of 0.54" diameter, corresponding to ~ 157 pc. Different emission lines and absorption bands are labeled. Blue 
labeled features are those used in this study. 
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the individual emission lines. We have used a nonlin- 
ear least-squares minimisation algorithm implemented in 
Image Data Language (hereafter, IDL) to perform the 
fitting. This method is comparable to other standard al- 
gorithms used for line fitting and implemented with other 
packages like IRAF. For each spaxel, local pseudocontin- 
uum levels were calculated for each spectral feature. Two 
bands were defined at each side of the emission lines and 
linear fits were applied between their average flux val- 
ues. In this way, our fitting algorithm weights each pixel 
with an error spectrum; it includes the Poisson noise, flux 
calibration uncertainties and errors in the local pseudo- 
continuum level. 

After subtracting the pseudocontinuum, the fitting al- 
gorithm returns fluxes, EWs and FWHMs among other 
parameters (like the Gaussian a and peak position in A) 
for all our lines. 

A similar treatment was used for the stellar absorp- 
tion features CO (2-0) (at 2.29 /xm) and Nal (at A2.206 - 
2.208 /im). Fluxes and EWs of the former were calculated 
by using the spectroscopic index definition of Kleinmann 
& Hall (1986). The Nal line features were measured by 
using the index definition from Forster Schreiber (2000). 
In the bottom panels of Fig.d] we present maps of the 
CO (2-0) feature; the weak Nal line is not suitable for 
mapping but for integrated region analysis only. Given 
the asymmetric shape of the CO (2-0) band and the dou- 
blet nature of the Nal absorption, we have not calcu- 
lated FWHMs for these features. We will leave that for 
Paper II, where the use of stellar libraries and spectral 
synthesis methods will be applied to derive the stellar 
kinematics. 

2.4. Flux, EW and FWHM Maps 

In Fig. [3] we present maps for the different phases of 
the interstellar medium (hereafter ISM) including the 
emission gas lines Br7 (at 2.17/im), the H2 transi- 
tion 1-0 S(l) (at 2.12 ^m), YielXlM [SiVI]A1.96 
and [FeII]A1.64 /im, including fluxes (none are corrected 
for internal extinction, see Sec. l3.2p . EWs and FWHMs 
(both in A). In a similar way. Fig. 2] shows the stellar 
components traced by H- and K-h&nA continua and the 
CO (2-0) band at 2.29 ^.m. (both flux and EW). All maps 
are centered close to the nuclear peak of the K-h&nd con- 
tinuum. At the distance of NGC 5135, our 8" x 8" FoV 
corresponds to the central 2.3 kpc x 2.3 kpc. The signal- 
to-noise (hereafter S/N) per A (per spaxel) is typically 
above 50 at the peaks of the emission, while in the dif- 
fuse regions it drops to about 20. The white areas on the 
maps (Fig.[3]and bottom panels of Fig. [4]) donnot present 
enough signal for reliable detections in individual spax- 
els. 

2.5. Selected Regions for study: Integrated spectra 

To better understand the processes taking place in the 
nucleus and neighbouring regions, we have chosen six 
separate regions for a more meticulous study. They have 
been labeled as A, B, C, D, E and F, and their loca- 
tions are indicated in the first row of Fig. [31 We refer the 
reader to Sec. 13. II for a detailed description of the rea- 
sons behind this selection. In Fig.[3l the circles' diame- 
ters correspond to the seeing of the observations within 
which we integrate our spectra to improve our S/N ratio. 



Such a S /N enhancement makes it possible to run certain 
tests and models with weak emission features (e.g. Nal, 
2-1 S(l)); this would be unfeasible on a spaxel- by-spaxel 
basis. At the distance of NGC 5135, we are integrating 
within aperture diameters of 157 pc in _ff-band and 180 pc 
in iiT-band. 

To extract the integrated spectra, we have used the IDL 
routine APER, available in the IDL Astronomy Library. 
This task allowed us to obtain fluxes from circular aper- 
tures by polynomial interpolation of the squared pixels. 
Aperture uncertainties were also provided by this pro- 
gram by using a close neighbours technic. As we did for 
the data-cube individual spaxels, we run our Gaussian 
fitting program on these integrated spectra, obtaining 
fluxes, EWs, FWHMs and other parameters, this time 
including the aperture error as an extra source of un- 
certainty. In Table[7| we present fluxes, EWs, FWHMs 
and their 1 a errors for each of our spectral features cor- 
responding to our six selected regions. In general, the 
flux uncertainties were between 10-20%, although weak 
features (e.g. 2-1 S(3)) have larger uncertainties. Typi- 
cally, the errors were dominated by Poisson uncertainties. 
However, for lines in the blue extreme of if-band (Br(5, 
1-0 S(3), [SiVI]) the flux calibration and pseudocontin- 
uum level determination are more uncertain given the 
variability and lower transmission of this spectral region. 

Finally, we performed internal flux extinction correc- 
tions for our measurements in regions A-F. Mean extinc- 
tion values were applied for the corresponding aperture 
spaxels derived from the extinction map (see Fig. [6] and 
Sec. 13. 21 for details). We present the corrected fluxes for 
all our lines in Table[7]with the mean Ay extinction ap- 
plied for each region. It is clear from the table that some 
of the uncertainties from the extinction corrections are 
relatively large compared to the uncorrected flux errors. 
In the extreme cases, the original <10% flux uncertainty 
of a strong line becomes an about 30% error in the cor- 
rected measurement. Therefore, we decided to use fully 
or non-corrected fluxes according to these criteria: Cor- 
rected fluxes were used when luminosities and other ab- 
solute flux dependant quantities (like hydrogen column 
densities, Nh) are calculated. Also, for the flux ratio 
[FeII]/Br7 we make use of fully-corrected fluxes for both 
lines. Non-corrected fluxes were only used for flux-ratio 
plots between lines in the TT-band. 

The reason for the criteria described above is to min- 
imise the uncertainties whenever possible. For the flux- 
ratio criteria, we found that the extinction correction be- 
tween the bluest and reddest lines in JT-band differs by 
6-8% depending on the A-F regions. In a flux ratio, this 
uncertainty is negligible compared to the individual flux 
errors in both lines. Therefore, while we do not signifi- 
cantly improve our measured ratio, we increase our error 
bars by about a factor of 3. When [Fell] A 1. 64 /^m is in- 
volved, however, the situation is different. In a similar 
analysis, we found that the extinction correction between 
[Fell] and a if-band line can easily differ by more than 
10%, even reaching 20% differences in the worst cases. 
These corrections are comparable to or larger than the 
errors of the non-corrected ffux ratios. Therefore, when 
[Fell] flux is involved, internal extinction must be taken 
into account. 

3. RESULTS AND DISCUSSION 
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Fig. 3. — Gas emission line maps of NGC 5135. White regions are those with unreliable detections for individual spaxels. 
From left to right: Br7, H2 (l-OS(l) transition at 2.12 /xm), HelA2.06^m [SiVI]A1.96 ^im and [FeII]A1.64 /im. From top to 
bottom: fluxes (non-corrected for internal extinction), EW and FWHM. The angular size of the spatial pixels corresponds to 
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In this section we present the results and analysis of the 
spatially resolved spectroscopy of NGC 5135. The many 
spectral features detected provide a variety of tracers of 
different physical phenomena. The most frequently used 
emission lines in the present work include Br7 (2.16 /im, 
recent star formation tracer), the H2 transition 1-0 S(l) 
(2.12 yum, warm molecular gas), [SiVI] (1.96 /xm, an AGN 
tracer) and [Fell] (1.64 /im, a supernova remnant tracer). 
Also, other features like the CO band at 2.29 /im, the Nal 
doublet at A2.206 - 2.208 ^m, HeIA2.06 ^m, Bt5 and the 
H2 transitions 1-0 S(0), 1-0 S(2), 1-0 S(3), 2-1 S(l) and 
2-1 S(3) were used during the analysis. 

3.1. Central Gas and Stellar Structure 

It is evident from the Figs.[3] and |4] flux maps that they 
are tracing different structures within the central 2.3 kpc 
of NGC 5135. In general, the number and spatial location 
of local flux peaks are totally different from one tracer to 
another. Only the nuclear region of this galaxy presents 
flux peaks in all maps. The EW maps (Fig. [31 second row 
and Fig.m bottom right panel) also differ substantially 
from tracer to tracer, showing local peaks usually match- 
ing those in flux. The only exceptions correspond to the 
nuclear region, where the AGN increases the continuum 
level (i. e. decreasing the EW values), and to the NE re- 
gion of the H2 map, where a EW peak does not have a 
flux peak counterpart. The FWHM maps (Fig. [31 third 
row) are rather smooth compared to the previous two. 
Only the FWHM ([Fell]) map presents clear structure, 
with two local peaks (north and south) separated by a 
third one in the nucleus. More detailed descriptions of 
some maps will be presented in the following pages, while 
some other enlightening comparisons between them will 
be used in different sections of this paper. 

As mentioned in Sec. l2.5[ the chosen regions A to F 
will help us to understand the different processes affect- 
ing the gas and stars in this galaxy. However, they will 
also be useful to label some of the structural features 
observed in Fig.[3l flux maps. Therefore, a brief descrip- 
tion of these regions follows. Region A corresponds to 
the nuclear peak in the if-band continuum. This region 
includes the AGN nucleus; all emission lines measured 
and the continuum have a peak at this location. Re- 
gions B, C and D correspond to bright peaks on the 
Br7 flux map, where probably very recent star forma- 
tion has taken place. Region E corresponds to the most 
prominent peak on the [Fell] flux map, presumably indi- 
cating the strong effects of shocks from supernova (SNe 
from now on). Also, this region matches the position of 
the secondary peak observed in the 1-0 S(l) transition 
of H2. We remark that the peaks of regions D and E, 
despite being close (0.56"), are not coincident. In terms 
of their integrated spectra, D and E share only a few 
common flux spaxels (equivalent to 1-2 spaxel's area) in 
the iT-band (0.62" aperture), while their _ff-band spec- 
tra are independent because of the smaller aperture used 
(0.54"). Any PSF-wing contamination effects were found 
to be negligible within these two apertures. Finally, re- 
gion F corresponds to a local peak on both the [Fell] flux 
and FWHM maps. It is also coincident with an extended 
structure in [SiVI]. In Fig.[l]and[2]wc show the extracted 
K and _ff-band spectra from these six regions. 

3.1.1. Coronal Gas 



The [SiVI]A1.96/.tm and [CaVIII]A2.32/.tm coronal hues 
are observed in the data, the former being our main AGN 
tracer (Fig. [31 fourth column; e.g. Maiolino et al. 2000, 
Rodriguez- Ardila et al. 2006a). The [CaVIII] hue is too 
faint and too close to CO (3-1) to be measured reli- 
ably. The ionization potential required to produce the 
[SiVI] hne (167 eV) is usually associated to Seyfert activ- 
ity where the gas is excited just outside the broad line 
regions of AGNs. 

As we can see in the Fig. [31 flux map, the [SiVI] line 
traces the galaxy nucleus and it also presents a weaker 
"plume" to the NE. This particular region, in terms of 
spatial scales, is at least 4x larger (w 600 pc in this 
case) than previous reports on different Seyfert galaxy 
samples (e.g. Prieto et al. 2005, Rodriguez- Ardila et al. 
2006a) where [SiVI] structures have at the most w 150 pc 
across. This finding, together with evidence for a simi- 
lar structure to the SW (see Sec. 13. 3. 2]) . suggest that we 
have detected, for the first time, the presence of ionizing 
cones in NGC 5135. According to Rodriguez- Ardila et 
al. (2006a), the morphology of [SiVI] and other coronal 
gas is preferably aligned with the direction of the tradi- 
tional lower-ionization cones (i.e. traced by [OIII]) seen 
in Seyfert galaxies. 

3.1.2. Ionized Gas 

The Br7 and Hel maps (Fig.[3l first and third columns) 
show emission dominated by three extra-nuclear sources 
at « 2" from the galaxy center, corresponding to re- 
gions B, C and D. As other recombination lines, the Br7 
feature is usually associated to recent star formation, 
where photoionization from massive young stars is tak- 
ing place. The higher ionization energy of Hel (24.6 eV) 
compared to hydrogen (13.6 eV) makes it a good tracer 
of the most massive 0-B stars and therefore, in princi- 
ple, of the youngest populations in starforming regions. 
An unambiguous interpretation of the Hel line, however, 
is rather controversial. Different authors have claimed 
that this emission depends on the nebular density, tem- 
perature and dust content, apart from the He/H rela- 
tive abundance and relative ionization fractions (Shields 
1993, Doherty et al. 1995, Lumsden et al. 2001). There- 
fore, without detailed photoionization models, certain 
frequently used tests like the Hel/Br7 ratio are very un- 
certain when analysing the hottest star temperature and 
tracing this population in Hii regions. This is why we 
will use Hel only as a first order spatial identification of 
the youngest stellar populations. The fluxes, EWs and 
FWHMs of this feature are included in Table[71 for com- 
pleteness. 

By comparing the Br7 and Hcl flux maps, it is evident 
that both are tracing the same structures, particularly 
the three extranuclear peaks mentioned above (i?, C and 
D) . It is worth noticing, however, that the strongest peak 
of the Hel-emission corresponds to the central region A, 
as it is also the case of Br7 after extinction correction 
(see Table[7]). These probably reflect that processes as- 
sociated to the AGN (e.g. shocks, photoionization. X-ray 
emission) contribute to the observed fluxes. However, 
given the size of our aperture (180 pc), we cannot dis- 
card the contribution of very recent starburst to further 
ionize the AGN vicinity. Gonzalez-Delgado et al. (1998) 
have shown for a small sample of local Seyfert 2 galax- 
ies that, even in UV/optical wavelengths (highly affected 
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by extinction), starforming knots are clearly identified at 
distances well below 100 pc from the nucleus (see also Cid 
Fernandes et al. 2004 and Davies et al. 2007). 

^From higher resolution (HST/NICMOS) Paa maps 
(Alonso-Herrero et al. 2006a) we know that at least re- 
gions B and C have a much more complex structure 
than observed here, and are composed of a few individ- 
ual knots of star formation. At least part of the impor- 
tant diffuse emission component observed in the Br7 map 
presents some structure in the Paa counterpart, partic- 
ularly within the "triangle" formed by ^, C and D. 

3.1.3. Partially Ionized Gas 

Different astronomical objects are capable of emit- 
ting [Fell] Al. 64 fim, including H ii regions, SNe remnants 
(hereafter SNRs) and AGNs, meaning that this line can 
be excited by different processes (see Sec. l3.3l for details). 

Strong [FeII]A1.64/im emission is produced by free 
electron collisions in hot gas (lO'^-lO''' K). Its low ioniza- 
tion energy (16.2 eV), however, indicates that this line is 
strong only in partially ionized regions; otherwise the Fe 
would be in higher ionization states. The SNR shocks are 
a suitable mechanism to produce partially ionized zones 
and therefore the [Fell] flux has been frequently used as 
a SNe rate indicator (e.g. Colina 1993, Alonso-Herrero 
et al. 2003, Labrie & Pritchet 2006, see Sec. 13. 61 for more 
references; see also Vanzi et al. 1996 for a different view). 
It is possible, however, that X-ray photoionization (e.g. 
from AGN or binary stars) produces a similar effect to 
that of shocks in the surrounding gas (Maloney et al. 
1996). 

Our [Fell] flux map is the most complex in structure 
(Fig-El fifth column), where two strong and three weaker 
extranuclear peaks are observed. The most prominent 
one corresponds to our E region and it seems to be 
a hot spot for SNe explosions given its high FWHM 
(« 513kms~^). On the other hand, our region F is also 
associated to a high FWHM zone but without a partic- 
ularly large [Fell] emission in comparison to the other 
extranuclear peaks. This region seems to match with 
the [SiVI] "plume" mentioned above. Developing more 
detailed kinematical arguments is beyond the scope of 
this paper and we will leave such analysis for Paper II. 
In Sec. l3.3.11 we will contrast the Fell with Br7 fluxes to 
shed some light on the dominant mechanisms responsible 
for heating the gas in the different regions. 

3.1.4. Warm Molecular Gas 

The observed H2 features in the spectra (Fig.[T]) are 
formed by roto-vibrational transitions of the molecules. 
Different mechanisms have been proposed as external 
energy sources including UV-fluorescence, thermal col- 
lisional excitation by SNe fast shocks and X-rays. 

Our strongest warm molecular gas tracer, 1-0 S(l), 
shows a very different structure than previous indica- 
tors (Fig. [HI second column). Two intensity peaks domi- 
nate the emission: the strongest at the nucleus and the 
weaker peaking in region E, matching the main peak in 
[Fell] flux. Both peaks are broad (the later, spatially 
resolved) , but the one in E is the dynamically hotter as 
shown in the FWHM map. The remaining H2 flux has a 
diffuse appearance at our spatial resolution. 




Fig. 4. — Stellar emission traced by diflerent indicators. 
White regions are those with unreliable measurements of CO 
in individual spaxels. Top panels: SINFONI K- and _ff-band 
continuum flux maps, the former indicating selected regions 
X, y and Z. Bottom: Flux of CO (2-0) band at 2.29 ^im (left 
panel) and its EW (right panel). The angular size of the 
spatial pixels corresponds to 0.125 x 0.125 arcsec^. 

In young starburst galaxies, massive (< 40 Mq) red gi- 
ants and supergiants with ages between 10-100 Myr dom- 
inate the integrated NIR continuum emission. As the 
population evolves with time, less massive giants start 
dominating the emission. Few Gyr later, the IR and 
bolometric luminosities are dominated by the low mass 
giants near the tip of the Red Giant Branch (Renzini & 
Buzzoni 1986, Chiosi et al. 1986, Origha & Oliva 2000). 
Spectral features such as the CO bands at 2.29 /xm and 
the Nal doublet (at 2.206 and 2.208 fj,m) are typical of 
K and later stellar types and also trace the luminous red 
giant and supergiant populations. The bottom panels of 
Fig.Hshow the flux and EW of the COA2.29^m band. 
White regions are those with unreliable detections of CO 
in individual spaxels. 

In Fig.m top panels, the H- and _K'-band continuum 
maps show a different spatial distribution than all pre- 
vious emission line maps in Fig.[3l We observed three 
extra-nuclear peaks which donnot correspond to any of 
our regions. However, from higher spatial resolution 
J-band continuum data (HST/NICMOS, Alonso-Herrero 
et al. 2006a) we know that these regions have a much 
more complex spatial structure with a series of knots. 
In any case, the strong diffuse emission component ob- 
served in their J-band continuum is also present in our 
SINFONI NIR data. 

In general, from the four panels of Fig. [5] we see that 
all these indicators are essentially tracing the same struc- 
tures, particularly in the extra-nuclear zone. The only 
exception occurs in the nuclear A region, where the 
strong NIR-continuum emission does not present a com- 
parable counterpart in COA2.29/im flux map. 

In Sec. l3.4l we will go further on characterising the stel- 
lar population and its spatial structure by using these 
tracers. We will also estimate its age in different regions 
by using simple stellar population models (SSP models 
from now on). 



3.1.5. Giant and Supergiant Stellar Component 



3.1.6. Overall Picture 
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Fig. 5. — Three contour maps overplotted: K-hand contin- 
uum (blue), [Fell] (black) and Br7 (red). In the same colors, 
local peaks of each map are marked by circles of approxi- 
mately the same size than our extraction apertures. 

Summarising, all the structural components described 
in this section clearly show different spatial distributions 
through their associated spectral feature maps. We wish 
to highlight that, apart from region A and region E (the 
[Fell] -main-peak / H2-secondary peak pair), many of the 
other line peaks are not spatially coincident with each 
other. We show this in Fig.[5l Three contour maps 
are overplotted including Br7, [Fell] and it'-band con- 
tinuum, where the circles represent the position of the 
peaks for each map, having approximately the diameter 
of our extraction apertures. According to StarBurst99 
models (hereafter SB99, Leitherer et al. 1999), these 
spectral features and iT-band emission might be trac- 
ing stellar populations on different evolutionary stages, 
pointing towards ages of « 6Myr, ^40Myr (assuming 
a SNR origin of [FeII]-emission) and ^200 Myr, respec- 
tively. In Secs. l3.3.T] and 13.41 we present a more de- 
tailed analysis of the stellar populations by using Br7, 
the COA2.29 urn ba nd, NaIA2.21^m and NIR-continuum. 
See also Secs. l3.3.T] and 13.61 for some discussion concern- 
ing the origin of [FeII]-emission in this galaxy. 

3.2. Internal Extinction Map 

Large amounts of dust and gas are usually found 
in local (U)LIRGs, hiding an important fraction of 
their star formation and AGN activity (Sanders et al. 
1991, Sanders & Mirabel 1996, Solomon et al. 1997). 
Dust grains are mainly responsible for the absorption 
of UV/visual photons and re-emitting them as sub- 
millimeter and FIR radiation. In the NIR range, dust 
also absorbs and scatters radiation, but its effects are 
much smaller than at shorter wavelengths. 

For NGC 5135, the effects of internal extinction be- 
come apparent in Gonzalez Delgado et al. (1998) work. 
By comparing their optical (606W at H5r/WFPC2, 
their Fig. 1) and UV-continuum maps {HST /FOG, their 



Figs. 5 and 6b) we found that our nuclear region A (cen- 
tered at +1.7"N, -t-0.4"W in their coordinate system) 
almost disappear in the UV compared to visual, while 
starforming regions like E retain their morphologies in 
both wavelengths. This suggests that the AGN neighbor- 
hood is more obscured than other circumnuclear regions 
of this galaxy. 

We can estimate the internal gas extinction in 
NGC 5135 by using observations of different hydrogen 
recombination lines. We adopt our own Br7 and Br(5 
measurements for this purpose. The theoretical ratio be- 
tween these two lines (Br7/Br,5 = 1.52 at T = 10 000 K 
and Tie = 10'*cm~'^, Osterbrock 1989) was compared 
with the measured fluxes for each spaxel. Therefore the 
extinction in magnitudes {A\) for both lines can be com- 
bined in the form 



A 



Br'j 



A 



BrS 



-2.5 



, (FBr-i,0\ 1 f FBr5,0\ 



= -2.5 • log 



{FBrf/FBrs)o 



{Fsr-i I FBrS)T 



(1) 



where F\^o and F\_t are the observed and theoreti- 
cal fluxes for a line centered at A, respectively. Then, 
by interpolating the extinction law of Calzetti et al. 
(2000), we put eq.[T] in function of Ay for each spaxel 
{Abt^ = 0.096 and AbtS = 0.132 Ay). From this 
expression we obtain the optical gas extinction map pre- 
sented in Fig. [HI White regions represent those where 
none flux information is available or where unreliable ra- 
tios were estimated. We remark that these extinction 
values are lower-limits only. This is because the gas 
might be optically thick (contrary to what is implicit in 
the theoretical values used above), therefore hiding the 
deeper regions from us. In terms of individual spaxels, 
the 1 a uncertainties in this Ay-gas map varies from 10- 
30 % to up to 80 %. The former typically corresponds to 
high S/N zones, like regions A-F, while the later matches 
with peripheral areas and some inter-region zones where 
the S/N is low. 

As we can see from Fig. [51 the regions with higher Ay 
values are: 

(i) The Nuclear Region: The average visual extinction 
is about 10 mag (0.9 mag at 2.2 ^m), with a local peak 
reaching the maximum of our scale at about 17 mag. This 
result is consistent with the conclusions of Gonzalez Del- 
gado et al. (1998). At least qualitatively, it is also in 
agreement with findings of Levenson et al. (2004) us- 
ing Chandra X-ray data. They claim that the nuclear 
AGN engine must be obscured by a high hydrogen col- 
umn density {Njj > 10^* cm~^, having in mind they have 
used 4x larger apertures than ours). Although we do not 
find such high Nh values (see Sec. [3X2l Tableg]), our 1- 
3 X 10^^ cm~^ measurements in region A are by far the 
highest of our A-F apertures. 

(ii) A region south from the nucleus (Aa w [—0.5, 1.0], 
A6 « [—2.0,-1.0], approx.): It has a mean visual ex- 
tinction of 10 mag. Presenting some local peaks of about 
17 mag, this area covers the space between our regions C 
and E. No strong emission in Br7, BrJ and Hel flux maps 
is observed in this area, but only a weaker, smooth com- 
ponent at SINFONI's resolution. Despite the fact that 
both Brackett lines are clearly detected in this zone, we 
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Fig. 6. — Ay gas extinction map of NGC 5135 based on 
Br7/Br5 line ratio. The arrows sliow our six selected regions 
A-F and the circles approximately represent the apertures in 
it-band. White regions are those with unreliable Av mea- 
surements in individual spaxels. The Calzetti et al. (2000) 
extinction law has been assumed. 

cannot discard that some of the peaks might be partiaUy 
produced by flux uncertainties, particularly of the BrJ 
line and its complex local continuum determination. A 
higher continuum level estimation would decrease this 
line flux, therefore, artificially increasing the extinction 
(see eq.[T]). In any case, because this region does not 
match with our A-F zones, it will not affect our extinc- 
tion corrected fluxes (TableO and conclusions derived 
from them. 

The rest of the map presents a median visual extinction 
of 6 mag which translates into 0.55 mag at 2.2 /im. 

Alonso-Herrero et al. (2006a) presented an FllOW- 
F160W color map of NGC 5135 (using NICM0S-NIC2; 
similar to J-H color) which can be interpreted as a stel- 
lar extinction map to a first order. It is difficult, how- 
ever, to directly translate this map to Ay extinction to 
the stars, mainly because we have to assume an age for 
the stellar population (to get its theoretical J-H color) 
in a galaxy which may have a composite stellar popu- 
lation (see Secs. l3.3.T] and In terms of their struc- 
ture, a comparison between Fig. [6] and Alonso-Herrero et 
al. color map shows that they are remarkably different. 
While both maps seem to agree in having the larger ex- 
tinctions around the nucleus, the color map clearly traces 
the nuclear spiral arms while no sign of this is observed 
in our Ay map. We do not think that the higher spatial 
resolution of NICMOS data can totally explain such a 
radical mismatch by itself. Instead, the close wavelength 
proximity between our two Brackett emission lines might 
be the responsible of the apparent Ay map lack of struc- 
ture. The similar extinctions for both Brackett emission 
lines produce Ay values (from eq.[T]) very sensitive to 
small flux uncertainties. Therefore, the detailed struc- 
ture in the Ay-ga,s map gets easily lost within Br7 and 
Br^ errors. We remark, however, that for larger areas 



(sampled by few tenths of spaxels like regions A-F), the 
mean Ay values should be a reasonably good estimation 
of the local internal extinction. 

3.3. Ionization Mechanisms of Gas 

The presence of ionized gas is usually interpreted as a 
tracer of recent star formation, where UV radiation from 
massive 0-B stars is capable of keeping the surrounding 
gas in an ionized state. In the NIR, recombination lines 
like Paa and Br7 are frequently used as tracers of this 
activity. 

However, there are other mechanisms capable of ion- 
izing the ISM. Shocks produced in SNR can partially 
ionize the gas. In the NIR, SNRs can be traced by 
the [FeII]A1.26 /zm and [FeII]A1.64/im lines (e.g. Green- 
house et al. 1991, Alonso-Herrero et al. 2003). A strong 
[Fell] flux is not expected in normal H ii regions (Mouri, 
Kawara & Taniguchi 2000): their abundant fully- ionized 
hydrogen would produce Fe atoms in higher ionization 
states. Therefore, by using our Br7 and [FeII]A1.64/j,m 
fluxes we are tracing different ionization mechanisms 
with also different ionization efficiencies. 

In Fig. [7] we present the extinction corrected [FeII]/Br7 
flux map. This diagnostic allows to see which ionization 
mechanism is more dominant in different regions of the 
galaxy. An [FeII]/Br7 ratio ^ 1 is usually taken as char- 
acteristic of a starburst, while much larger values are 
associated with only partially ionized regions by shocks 
(e.g. Graham et al. 1987, Kawara et al. 1988, Mouri et 
al. 1990, Alonso-Herrero et al. 1997, 2001). 

3.3.1. Circumnuclear Starforming Regions 

Fig. [7] shows that only three regions (associated to B, 
C and D) seem to have a non-negligible contribution 
from starburst activity. In the same figure, we overplot 
in contours the EW(Br7) as an upper-limit age indica- 
tor. It is clear that the EW(Br7) peaks match very well 
those areas with low- [Fell] /Br7 ratio, pointing out that a 
young stellar population contributes to ionize the ISM in 
these areas. We have tested this idea by using the SB99 
models with instantaneous star formation, solar metal- 
licity and Kroupa initial mass function (hereafter IMF) . 
In Table[l] we present the predicted ages (upper limits) 
of the stellar populations for our selected regions, de- 
rived from our EW(Br7) measurements of Table[7l As 
we can see, B, C and L) are the three youngest areas 
with w 6Myr old, while E and F are slightly older. 
These results are consistent with previous findings with 
UV/visual data (Gonzalez Delgado et al. 1998). It is 
possible that some contamination from the nearby star- 
burst regions has been introduced within E and F zones, 
particularly for the former. ^ 

A plausible explanation for the local [FeII]/Br7 varia- 
tions might be suggested from Table[T]ages, where regions 
E and F are marginally older than the others. According 
to models, after about 6-7 Myr important changes start 
taking place in young stellar populations. The massive 

The age presented for region A is an upper limit for two reasons: 
as for the other regions, we are assuming that all Br7-emission 
comes from star formation; but also in region A the presence of 
the AGN increases the continuum level producing lower EW(Br7) 
values and therefore, larger ages than expected for a pure starburst 
area. 
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StarBurst 99 Ages from EW(Br7) 
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"•^^-0.03 
fi 51+0-04 



(") Ages with their itlcr errors. Ages should be considered as upper 
hmits only. See text for details and footnote 3 for region A. 
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Fig. 7. — Extinction corrected [FeII]/Br7 line ratio map. 
In contours we plot the EW(Br7) as an age indicator. The 
arrows show our six selected regions A-F and the circles ap- 
proximately represent the apertures in if-band. 

0-B stars responsible for most of the Br7-emission ex- 
plode as SNe, diminishing this flux feature in the inte- 
grated spectra. On the other hand, after 4 Myr since 
the star burst, the SNe rate has been increasing dramat- 
ically, enriching the ISM with Fe. Therefore, two dif- 
ferent processes radically change the [Fe/H] abundance 
and the Br7 flux: when total ionization from 0-B stars 
leaves place to partial ionization from SNR shocks, the 
Fe-enriched ISM increases the [Fell] flux, while the de- 
cline of 0-B stars decreases Br7. If a certain region has 
not reached yet this crucial stage on its stellar evolution 
(like B, Cand D) we would observe low [FeII]/Br7 ratios. 
On the other hand, if the stellar population has reached 
this point (like in E and F), much higher [FeII]/Br7 ra- 
tios would be observed. 

The high values of [FeII]/Br7 (> 20) in regions iJand F 
are, in fact, very similar to those observed in SNRs, e.g. 
- 30 in IC 443 (Graham et al. 1987) or 34 in RCW 103 
(Oliva et al. 1989). In the FWHM([FeII]) map of Fig.[3] 
we see that both sources are the only ones with very high 
values (or large velocity dispersions). This would be, in 
principle, in excellent agreement with a scenario where 
SNR shocks have perturbed the kinematics of these areas. 

We have compared all our [FeII]/Br7 measurements 
with predictions from the SB99 models. On its latest 
versions, SB99 does not give directly the [Fell] fluxes 



O Region A 

+ Region B 

* Region C 

A Region D 

□ Region E 

X Region F 



1.0 



0.5 



0.0 




6.76 



6.78 



6.80 6.82 
log(Time) [yr] 



6.84 



6.86 



Fig. 8. — [FeII]/Br7 vs. Age (upper limits from EW(Br7)) 
for our six regions. Based on SB99 models, different relations 
have been used to derive the [Fell] flux from the SNe rate: 
dashed line, Calzetti (1997); solid line, Alonso-Herrero et al. 
(2003). 

but the SNe rate. So we transform these values to [Fell] 
flux by using empirical relations from Calzetti (1997) and 
Alonso-Herrero et al. (2003). By using these expressions, 
we can further test (to first order) the SNR scenario. 
See Sec. 13. 61 for further discussion about SNe rate esti- 
mations. In Fig. [8] we present SB99 predictions with the 
two empirical relations mentioned above. We overplot 
our measured [FeII]/Br7 ratios for ^-Fwith their uncer- 
tainties. We have used the upper-limit ages of Table[T]to 
flx the abscissa of our data points. 

Only regions D and F show a clear agreement with 
the SNR models at those ages. Region D has an im- 
portant [Fell] emission contribution from the broad peak 
traced by region E. So, the agreement between this young 
starburst region and the models is not surprising. The 
good match with F only reinforces the link between this 
area and SNRs as mentioned before using kinematic ar- 
guments. The starburst B and C regions lay below model 
predictions, pointing towards pure stellar photodissocia- 
tion mechanisms as the responsible of the gas ionization. 

Certainly, the most interesting outlier is region i?, 
where the models predict a smaller [Fell] /Br7 ratio than 
observed. This implies that an additional ionization 
source must be invoked to explain the observations at 
the center of this resolved [Fell] peak (region D, taken as 
a peripheral sample of this peak, does not need such an 
extra ionization source). A plausible alternative comes 
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from X-ray photoionization. As mentioned in previous 
sections, Levenson et al. (2004) have published Chan- 
dra X-ray data for NGC5135. They have identified in 
their total band map (0.4-8 keV) a strong emission cor- 
responding to region E. According to these authors, the 
non-thermal component (power law and the bulk of the 
hard X-ray emission) is likely to be associated to X-ray 
binaries in the starburst, while the thermal one (most of 
the soft emission) is typical of a starburst galaxy. These 
mechanisms may explain the difference between models 
and iJdata in Fig. [51 where both, X-rays and SNR shocks 
are contributing to the observed [Fell]/Br7 ratio. 

3.3.2. The Ionization Cones and Central ACN 

The central region A of NGC 5135 deserves a separate 
section given its particular characteristics with respect 
to the other starforming regions. 

As mentioned in previous sections, the AGN activity 
has been traced by the [SiVl] coronal line. The [SiVll- 
luminosity for this region is L[Sivi] — 3-6 x 10^^ ergs" , 
in agreement with similar measurements in other local 
Seyfert2 galaxies, like Circinus, NGC 1068, NGC 3227 
and Mrk609 which span a range ~ 10'^^-10'*° ergs^^ (the 
two later objects are Seyfert 1.5 and 1.5-1.8, respectively; 
Rodriguez- Ardila et al. 2006a, Zuther et al. 2007). 

By using Br7-emission, we can make a row estimation 
of the AGN contribution to the gas photoionization. If 
we assumed that the (extinction corrected) Br7 flux in 
region A is entirely produced by the AGN, we can com- 
pare this value with the total Br7 flux in our FoV. Ac- 
cording to this, the AGN contributes in < 3% of the gas 
photoionization within the central 2.3 kpc of this galaxy. 
By doing a similar calculation, we also estimate that the 
AGN represents < 3 and < 5% of the NIR-luminosity in 
the H- and JT-bands, respectively. 

Concerning the ionization mechanisms, on the one 
hand the AGN embedded in this zone is likely to play an 
important role at local level by photoionizing the gas in 
the nuclear neighborhood (Marconi et al. 1994, Ferguson 
et al. 1997, Contini et al. 1998, see also Rodriguez- Ardila 
et al. 2006b and references therein). The [FeII]/Br7 = 5.7 
measured in this region is consistent with typical Seyfert 
galaxy values from a compilation of different object types 
shown in Alonso-Herrero et al. (1997). 

The "plume" shape structure in [Si VI] -emission (NE 
of the nucleus. Sec. 13.1.11 Fig. [21) probably traces the di- 
rection of a cone, pointing that the AGN is also capable 
to ionize the ISM at larger distances. The presence of 
shocks from outflowing winds would produce such ionized 
plasma and cone-features (Viegas-Aldrovandi & Contini 
1989). The [Si VI] -emission is particularly broad in A 
with a FWHM of « 28 A (or « 430 km s"^), strengthen- 
ing the idea of an important interaction between AGN 
outflows and the surrounding material. It is interesting 
to notice that the FWIIM([SiVI]) nuclear peak structure 
extends approximately along the semi-minor axis as ex- 
pected in outflow models of Heckman et al. (1990). This 
may suggest the presence of multiple components in the 
[SiVI] kinematics. The high FWIIM(Br7) in this area 
(« 22 A or « 300kms~^) is also consistent with this hy- 
pothesis. In Paper II we will find out if this scenario 
receives further support from our detailed kinematical 
study. The NE [SiVI] structure has a SW counterpart 
not detected on a spaxel-by-spaxel basis. By integrating 



few spaxels at 2" SW from the nucleus, we have detected 
clear [SiVI]-emission which is absent in other areas close 
to region A. This SW component is probably weaker be- 
cause of extinction and projection effects. 

On the other hand, in Levenson et al. (2004) work on 
X-rays, the authors have found the strongest emission 
peak in our region A, which may also contribute to the 
ionization in this area. They suggest, however, that the 
bulk of the thermal (coUisional) X-ray emission, even in 
the nucleus, is stellar^ not due to the AGN. In principle, 
this would be understood in terms of the larger apertures 
they use in A (about 4 x ours) within which a stronger 
stellar component, that we miss, might be introduced. 
However it is also true that compact starforming regions 
have been usually found in AGN's neighborhood at dis- 
tances of tens of pc (e.g. Gonzalez Delgado et al. 1998, 
Weaver 2001, Davies et al. 2007). Considering the signif- 
icant dilution from the AGN, the later point is reinforced 
for NGC 5135 by our findings of Nal and CO stellar ab- 
sorption in this zone. 

Therefore in region A and in the [SiVI] -cone the 
AGN photoionization and outflows seem to play an 
important role on ionizing the gas. However, a stellar 
contribution cannot be totally discarded, particularly in 
region A. Then, it is comprehensible why in Fig.[8l the 
starburst-based models cannot reproduce the observed 
[FcII]/Br7 ratio of this region, in addition to the age 
uncertainty in this zone. 

Finally, the remaining extended emission in Fig. [71 
has rather large [FeII]/Br7 ratios, typically between 3 
and 7. This suggests that the ISM gas is mainly ionized 
by other mechanisms (e.g. SNR shocks. X-rays) instead 
of photodissociation from star formation. 

3.4. Dating the Stellar Population: Further constrains 
from Absorption Lines 

In the previous section, we have studied a young stel- 
lar population (< 7Myr) with either an important frac- 
tion of 0-B stars or in a transition period with high SNe 
rates. We attempted to understand its importance on 
ionizing the ISM gas by both, the use of the [FeII]/Br7 
map structure, and the estimation of the stellar popula- 
tion age. Because of this, we now consider appropriate 
to introduce at this point the analysis of the stellar pop- 
ulation traced by the NIR-continuum and COA2.29/im 
absorption maps of Fig.[4l By doing this, we will address 
the open issues presented in Sec. 13.1.51 

In Fig.m we have seen that NIR-continuum and 
COA2.29/xm maps have a very similar structure every- 
where but in region A. The absence of a strong central 
peak in the CO fiux map can be explained by the influ- 
ence of the AGN in the stellar continuum. The strong 
X-ray/UV radiation from AGN is reprocessed by the sur- 
rounding dust grains and re-emitted in the IR (e.g. Clavel 
et al. 1992, Forbes et al. 1992). Therefore, the contin- 
uum fiux in H and JT-bands increases, diluting the CO 
and other stellar absorption features. The high extinc- 
tion we estimate for this area (see also Alonso-Herrero et 
al. 2003) and the finding of hot dust emission (^ 200 K) 
in the nucleus of NGC 5135 (Gemini/T-ReCS A^-band 
data, Alonso-Hererro et al. 2006b) are consistent with 
this idea. The high temperatures required for the dust 
to emit in JT-band, however (1200-1500 K; e.g. Clavel 
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et al. 1989, Glass 1992, Rodriguez- Ardila et al. 2005a, 
2006b), limit this effect to rather short distances from 
the AGN (< 100 pc; Marco & AUoin 1998, 2000). In fact, 
the observed dilution towards the center of the EW(CO) 
map (Fig.m bottom-right panel) occurs at similar dis- 
tances of about 0.4" (~ 120 pc). This apparently close 
agreement, however, is almost certainly driven by PSF 
effects. It is possible than hot dust emission is confined to 
even shorter distances from the active nucleus. In regions 
farther away from the nucleus (> 200 pc) the EW(CO) 
values are free of the AGN contamination and can be 
confidently used for further tests of the stellar popula- 
tion properties. 

As for A-F zones, we extracted integrated spectra for 
other three extra-nuclear regions labeled as X, Y and Z 
(see Fig.m top-left panel). They do not match with our 
A-F apertures but coincide with the three extra-nuclear 
peaks of the NIR-continuum (this is clearly shown in 
FiglS]) where the CO signal is higher. Our aim is to do a 
general comparison between the young and (potentially) 
old stellar populations for the entire galaxy, based on in- 
dividual regions where one or the other is dominant. In 
this way, we base our conclusions on higher S/N data. 

We estimate the stellar population ages in X-Zhy using 
the EW(COA2.29 fim) and the models SB99 and STARS 
(Sternberg 1998, Thornley et al. 2000, Davies et al. 2003). 
Both models include treatments for the thermally puls- 
ing AGE stars. However, without knowing parameters 
like the metallicity and characteristics of the star forma- 
tion history, the model predictions might vary from case 
to case. We have run these models using different IMFs 
(Kroupa, Salpeter) and metallicities (0.04, 0.02, 0.008) 
without finding substantial variations which may change 
our conclusions. Therefore we have assumed solar metal- 
licity, Kroupa IMF and the Genova isochrons for all our 
age predictions. Only the use of different star formation 
histories (instantaneous and continuous) has a relevant 
impact on this respect. We will discuss this particular 
point in the following paragraphs. 

Strictly speaking, STARS makes use of exponentially 
decaying star formation rates (hereafter SFRs). How- 
ever, by using a long characteristic time-scale (1000 Gyr) 
we recover a continuous SFR, while by using a short 
one we basically recover instantaneous burst results for 
ages > 1 Myr. We can directly compare this models to 
instantaneous and continuous SFR models from SB99 
(R.L. Davies, private communication). 

In Table[2] we present instantaneous-burst model re- 
sults for regions X-Z. Because SB99 and STARS models 
interpolate the EWs from spectral grids, their predictions 
become quite noisy at larger ages (> 100 Myr), particu- 
larly for STARS were a template star library is used (see 
FiglU R.L. Davies, private communication). Therefore, 
in cases of old stellar populations (> 100 Myr) a range of 
ages is presented for SB99, while lower limits are shown 
when STARS is used. 

As we can see in Figl9] (upper panel), instantaneous- 
burst models produce two possible stellar population 
ages consistent with our data (around 7 and > 200 Myr, 
respectively), while continuous-star- forming ones pro- 
duce a single range of ages around 7 Myr. Although a 
stellar population of 7 Myr would agree with our find- 
ings in 5-F zones using EW(Br7), the presence of a sec- 
ond, totally different population receives support from 
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Fig. 9.— SB99 and STARS model predictions for EW(CO) 
and EW(Nal). The two horizontal lines represent the range 
of our data for regions A-F and X-Z. 

two different sources. 

First of all, our CO-fiux (Fig. [4] bottom left panel) and 
Br7-emission maps (Fig. [3] top of first column) clearly 
trace different structures (peaks) in this galaxy, usu- 
ally separated by more than 200 pc in projection. On 
the other hand, the CO-flux map almost mimics the 
structure observed in the NIR-continuum, which in turn 
is frequently used as a tracer of more evolved stellar 
populations than those with high EW(Br7). Further 
support for this argument comes from Alonso-Herrero 
et al. (2006a) NICMOS/ffST data of this galaxy (see 
Alonso-Hcrrero et al. 2006b, Fig. 1 second row for sharper 
figures). The higher spatial resolution of their NIR- 
continuum and Paa-emission maps show the peaks ob- 
served with SINFONI actually have a complex struc- 
ture. Moreover, the NICMOS continuum and Paa peaks 
are not coincident spatially, demonstrating that, between 
the neighbouring aperture pairs A- 5, Y-C and Z-D^ we 
are picking-up different spatially resolved regions (even 
though some PSF-contamination effects between aper- 
ture pairs). In terms of structure, further support to spa- 
tial variations in the stellar population comes from the 
EW(CO) map (Fig. |4] bottom right panel). This figure 
shows a clear radial (projected) distinction in function 
of EW. While at distances of ^ 600 pc from the nucleus 
we find EW(CO) 7A, this value decreases to ~ SA 
at ^ 300 pc. Such spatial segregation supports the idea 
of a change in the stellar population properties (e.g. a 
variation in the relative contribution of red giants and 
supergiants to the integrated hght). 
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TABLE 2 

COA2.3/.tM, NAlA2.21/iM and Br7 EWs and Instantaneous burst model ages for regions X-Z 



Reg 


EW(CO) 

[A] 


[Myr] 


[Myr] 


^&''^SB99 

[Myr] 


[Myr] 


EW(Nal) 

[A] 


Agcl^'TMs 
[Myr] 


Agc2^';^;3 

[Myr] 


EW(Br7) 

[A] 


^g^SB99 

[Myr] 


X 
Y 
Z 


6.44 ± 0.4 
7.14 ± 0.3 
6.90 ± 0.3 


7 q7+u.ua 

' -0.03 
7 42+0.02 
' -^^-0.02 

'•^'^-0.02 


'•'-'^-0.08 
7 14+0.05 

7 -|n + 0.05 

' •^'J-o.io 


213-453 
205-211 
205-445 


> 200 

> 200 

> 200 


1.76 ± 0.08 
1.99 ± 0.07 
1.88 ± 0.08 


o OO-I-0.04 
*-^*-0.06 

».^y_0.06 


> 200 

> 200 

> 200 


6.01 ± 0.64 
7.15 ± 0.53 
6.64 ± 0.57 


6.70lH^ 



Secondly, the EW of different spectral features also 
point towards different stellar populations in X-Z with 
respect to B-D zones. On the one hand, the EW(CO) of 
the former are larger than those for the later as we can 
see in Tables[5]and[71 respectively. The quoted uncertain- 
ties show that these differences are significant at a 4-5 (t 
level. We observe the EW differences are mainly driven 
by CO-absorption variations rather than continuum-level 
changes. This suggests the existence of a different stel- 
lar population in X-Z with respect to B-D^ the former 
regions having a more important red giant/supergiant 
contribution to the integrated spectra (e. g. Origlia & 
Oliva 2000). On the other hand, as we can see in Tables[2] 
and [71 the EW(Br7) is substantially lower in regions X- 
Z than in B-D zones, also pointing towards two different 
populations. Putting together all the observational ev- 
idence, the NIR-continuum and CO-absorption (red gi- 
ant/supergiant stars) and Br7-emission (0-B stars) sug- 
gest the existence of (at least) two stellar populations in 
this galaxy. 

However, as mentioned above, the stellar population 
model predictions are not that clear about this point, de- 
pending on the star formation history used (see Figl9]). 
If we assume continuous SFR models, the EW(CO) for 
all our apertures are consistent with a 7 IVIyr old pop- 
ulation. By using the EW(Br7) with the same mod- 
els, however, we get different results. The B-D regions 
are now consistent with a 10 IVIyr old stellar population, 
while regions X-Z are better constrained by a ~ GOlVIyr 
old one. Clearly, the use of continuous models with the 
measured EWs produce inconsistent results between the 
two tracers. 

Instead, if we use instantaneous SFR models, the ob- 
served EW(CO) are consistent with two populations of 
7-8 and > 200 Myr old, respectively (see Table[2] and [3] 
for regions X-Z a.nd B-F, respectively). This degeneracy 
is broken by using EW(Br7). Even though regions X-Z 
have much smaller EW(Br7) values than B-D, all these 
zones are consistent with a 6-7]VIyr stellar population 
according to models (Tables^] and [3]). In this case, the 
instantaneous model predictions are consistent between 
our different EW measurements. 

We have also compared the EW(NalA2.21/^m) mea- 
surements with STARS model predictions. All previous 
warnings and discussion about the EW(CO) age estima- 
tion are also applicable to this stellar absorption feature. 
As we can see in FiglS] bottom panel and in Tables^] 
and [3l the younger-age predictions for an instantaneous 
model are somehow larger (~ SlVIyr) than those with 
EW(CO), but still consistent with a young population. 

In summary, while observational evidence points to- 
wards two different stellar populations, the model pre- 
dictions cannot clearly discriminate between them. Our 
measurements are consistent with temporally close, in- 



stantaneous bursts occurred 6-8 Myr ago. Strictly speak- 
ing, regions X-Z are always older than their neighbors 
B-D (independently of the tracer used) by few hundred 
thousand years. Even though it is worth mentioning, 
such small age differences might be in the limit of cur- 
rent model capabilities for discriminating between two 
stellar populations. Certainly, we cannot discard the pos- 
sibility of two populations of similar ages: a young one 
(~ 7 Myr) dominated by 0-B stars and an intermediate- 
age one (~ 10 Myr) with massive stars which already 
started their evolution to red supergiants. In such partic- 
ular case, the presence of local metallicity variations be- 
tween regions might become an issue given the tight cor- 
relation between the CO-index and metallicity (Origlia 
& Oliva 2000). In any case, given the above discrepan- 
cies, we cannot be conclusive about the existence of one 
or more stellar populations and thus, further data and 
tests are necessary to shed light about this issue. 

3.5. H2 Excitation Mechanisms 

Different mechanisms have been proposed to be re- 
sponsible of the H2 excitation and literature is rich 
in examples where more than one process is shaping 
the observed emission spectra (e.g. Draine & Woods 
1990; Mouri 1994; Davies, Sugai & Ward 1997; Quillen 
et al. 1999; Burston, Ward & Davies 2001; Davies et 
al. 2005; Rodrfguez-Ardila et al. 2005b; Zuther et al. 
2007). Between these mechanisms, the most commonly 
discussed are UV- fluorescence (photons, non-thermal), 
shock fronts (collisions, thermal) and X-ray excitation 
(collisions through secondary electrons, thermal). 

In this section, we will focus on trying to understand 
the observed ll2-emission by using different diagnostic di- 
agrams. We will consider both non-thermal and thermal 
excitation mechanisms and models in the analysis. 

3.5.1. H2 Line Ratios: A Diagnostic Diagram 

As a first step, we have calculated different H2-emission 
line ratios for our six regions in NGC 5135. As mentioned 
in Sec. l2.5i we did not use fluxes corrected for extinction. 
Different models and data suggest that the excitation 
mechanisms mentioned above produce characteristic line 
ratios between different roto-vibrational transitions (see 
Fig. [To)) . In this sense, the many lines presented in our 
K-hand data allowed us to use a line-ratio diagnostic 
diagram for our analysis. We maximise the difference 
between thermal and non-thermal processes by trying 
to remove the ortho/para ratio effect (ratio equal to 3 
for a thermal distribution and < 3 for a non-thermal 
one). We adopt the intensity ratios 2-1 S(l)/l-0 S(l) 
and 1-0S(3)/1-0S(1) to achieve this goal. Given that 
the 2-1 S(l), 1-0 S(3) and 1-0 S(l) lines occur in ortho- 
molecules only, the non-thermal values of their ratios are 
supposed to be constant (Mouri 1994). 
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TABLE 3 

COA2.3/^M AND NaIA2.21^m Instantaneous burst model ages for regions B-F 



Reg 


Agel^^gc, J^^L Ago2^^^'gg Agc23V,, 
[Myr] [Myr] [Myr] [Myr] 


[Myr] [Myr] 


C 
E 

pa 


7-26+0 02 6.72+^;^^ 534-592 > 811 
7.33to E;2 6.9llo:°g 217-509 > 200 
'i'-SOlo.oa 6.82l!5:o7 233-552 > 811 
7.4llg g^ 7.12l!5:o7 200-443 > 200 
7.22^0 02 6.6llo;°7 581-698 > 811 


7.89+0 05 ^ 200 
8.05to g| > 200 
8.1ll0;07 > 200 

8.39t°;o7 > 200 
S.Olto og > 200 



" The apparent inconsistency between different model predictions for regions B, D and F oldest ages is mainly driven by the noisy nature 
of these models for ages > 100 Myr. 



In Fig. [To] we present the line-ratio diagnostic diagram 
for our six regions compared with different models and 
data from literature. The first obvious conclusion from 
this diagram is that none of our six regions seem to be 
dominated by non-thermal (fluorescence) excitation. In- 
stead, they lie closer to thermal processes like X-ray heat- 
ing and SNR shocks, the most likely dominant mecha- 
nisms responsible of H2 excitation in the six zones. Re- 
gions B and D, however, clearly depart from the pure 
thermal emission curve suggesting some UV-fluorescence 
contribution to the observed lines. This is reasonable 
in principle given that B and D correspond to recent 
star forming regions where UV-photons are abundant. 
The remaining four regions C, E and F) are consis- 
tent with the purely thermal emission curve (within the 
uncertainties) with gas temperatures between 2000 and 
3000 K. The overall lack of dependence on UV-processes 
(i. e. recent star formation) agrees with the comparison 
between our Br7 and H2 flux maps (Fig. [J]). The clear 
structural differences between them implies that the pro- 
cesses exciting the H2 gas are not those who produce Br7 
(i.e. recent star formation). 

By comparing with models, it is interesting to notice 
that region A, the AGN nucleus, matches very well with 
the X-ray models of Draine & Woods (1990). As men- 
tioned in previous sections, Levenson et al. (2004) have 
found a strong X-ray source in this area. * In their dis- 
cussion, Draine & Woods (1990, 1991) argued that there 
is a range of values of X-ray energy per volume for which 
the gas gets sufficiently hot to vibrationally excite the H2 
(producing the observed II2 lines with high enough effi- 
ciency to be observed), but not hot enough to dissociate 
it. They proposed the X-ray emission from SNRs (and 
not their shocks) as one plausible stellar energy source. 
This scenario implies other constraints like high gas den- 
sities {uh > lO^cm"^, attainable in the dense circum- 
nuclear regions around AGNs) and high SNe rates (e.g. 
« 1.5yr-i within the central 4kpc of NGG6240, but 
with high uncertainties). The later constraint seems to 
be achieved by NGC5135, with a SNe rate of ~ 0.4- 
0.6 yr~^ within the central 2.3 kpc (see Sec. l3.6l for fur- 
ther details). As an extra possible piece of evidence, the 
Draine & Woods' models predict a decrement in the 2- 
1S(3) transition flux (w = 2, J = 5 level), which we 

^ In Levenson et al. (2004) the aperture's diameters are about 
4 X those of our selected regions. We have repeated all the tests 
of this section by using Levenson et al. (2004) aperture size and 
no significant changes have been found. The relation between data 
and all model predictions reminds the same, without affecting our 
discussion. 
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Fig. 10.— H2 diagnostic diagram using the 2-1 S(l) /l-O S(l) 
and 1-0 S(3)/l-0 S(l) intensity line ratios. Diamonds repre- 
sent our own data for the six selected regions of NGC 5135. 
The X symbols represent thermally (Sternberg & Dalgarno 
1989) and non-thermally (Black & van Dishoeck 1987) UV- 
exited models. Squares correspond to X-ray heating models 
of Draine & Woods (1990). The asterisk is a shock model 
from Brand et al. (1989). The triangle is SNR data from 
Oliva, Moorwood & Danziger (1990). The curve represents 
thermal emission between and 3000 K (eqs. [2] and [3] see 
references in the text). 

actually observed in Fig.[TT] (within large uncertainties). 

We also highlight in Fig.fTU] that region E is consistent 
with a thermal process but far from the SNR data and 
shock model. As we saw in Sec. 13. 31 SNR shocks were an 
important candidates for Il2-excitation given the high 
[Fell] fluxes and FWHMs observed in this area. Could 
region i? be in a similar scenario to the one described 
above for Al Certainly, the X-ray models are not as close 
to E as to region ^4, however we know from Levenson et 
al. (2004) work that the second X-ray peak comes from 
region E and it has a star formation origin. In principle, 
the same arguments presented before for region A might 
be applied to E in the sense that an important fraction 
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of H2 excitation might be produced by X-rays of stellar 
origin. 

Finally, region F is also consistent with thermal pro- 
cesses. However, the uncertainties make it difficult to 
distinguish between X-ray and SNR shock scenarios. Be- 
cause the later has a dominant role in ionizing the local 
ISM gas fSec. l3.3.r|) . SNR shocks become a likely candi- 
date to dominate the H2 local excitation in this area. 

3.5.2. H2 population diagram 

Because we are able to measure H2 -emission lines from 
different vibrational transitions (1-0 and 2-1), we adopt 
another approach to characterise the observed excita- 
tion mechanisms. We have just seen in Sec. 13.5. 1] that 
the dominant processes responsible for H2 excitation are 
thermal (coUisional). Therefore, if we assume pure ther- 
mal excitation, the different roto-vibrational levels are 
populated as described in the Boltzmann equation. This 
allowed us to represent the different transitions as a func- 
tion of their population (probability) densities versus the 
energy of the upper level, in what is called a population 
diagram. For a thermally excited gas, all the transition 
values lie on a straight line and the corresponding slope 
is inversely proportional to the H2 gas temperature. 

Given that molecular gas clouds contain a large num- 
ber of molecules, the population densities are essentially 
proportional to the observed column densities. So, the 
probability ratio between two different states can be ex- 
pressed with the Boltzmann equation as 

Pa ^ N{v,J)a ^ gg ■ 6-^°/^^ 

Pfc N{v, J)b gb ■ e-E^I^^ ^ ^ 

where Pa and Pf, are the probabilities of finding an H- 
molecule in any of the degenerate states (ga, Qh) with 
energy Ea, P&. N(v, J)a and N{v, J)h are the observed 
column densities at the corresponding states, v and J 
are the vibrational and rotational quantum numbers, re- 
spectively, k is the Boltzmann constant and T is the gas 
temperature. The column densities can be derived with 
the formula: 

N{v,J)^ / , ,n x-Ax-^ (3) 

A{vJ,v'J') he ilaper 

where / is the measured flux for the corresponding emis- 
sion line, A{v J, v' J') is the transition probability from 
the {v J) to the {v' J') quantum state (taken from Wol- 
niewicz et al. 1998), A is the rest frame line wavelength, 
h is the Planck constant, c is the speed of light and f^aper 
is the squared aperture diameter of 0."62. The derived 
column densities for each H2 indicator and region are 
presented in TableSl 

By taking the logarithm in eq.[2] we end up with a 
generic expression for N{v, J) of the form 

f N{v,J)/gj \ -E{v,J)/k , ^ 

^" U(l,3)/g3 J = f + 

where we have chosen the transition {v, J) — (1,3) (cor- 
responding to the 1-0 S(l) line) to normalise the derived 
column densities. The constant is independent of the 
transition. Because the ratios N{v, J) /Nil, 3) are equiv- 
alent to flux ratios, we have used the none extinction 
corrected fluxes for the population diagram as argued in 
Sec. 12. 51 In any case, we present in Table|4] the column 
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Energy [K] 

Fig. 11. — H2 level population diagram relative to l-OS(l) 
for the six regions of NGC 5135. The data points are slightly 
spread in the X-direction simply for illustrative purposes. The 
straight line is the linear fit to the v — 1-0 levels. Its slope is 
consistent with an isothermal population at 1930 ± 114 K. 

densities derived from each H2 line with and without ex- 
tinction correction. 

In general, our corrected column densities, Nucon 
range between ~ 10^^-10^^ cm~^ for all indicators and 
regions, with slightly lower densities derived from the 
V — 2-1 transitions. The larger Nu.corr values coming 
from region A (1-3 x 10^^ cm~^) are not as large as those 
derived by Levenson et al. (2004) (> lO^^cm^^, from X- 
ray observations). Gonzalez Delgado et al. (1998, using 
UV data) and Levenson et al. (2004) have also reported 
column densities for E (9.9 x 10^° and 1.1 x 10^^ cm~^, 
respectively but within larger apertures than ours) which 
are consistent with many of our estimations in this area. 

In Fig.[Tl]the II2 population diagram for our six regions 
is presented, showing the column densities normalised by 
the (1,3) transition versus the energy of the upper level 
(in Kelvin). We have used six different emission lines 
including four v — 1-0 transitions (S(0), S(l), S(2) and 
S(3)) and two v = 2-1 transitions (S(l) and S(3)). The 
different symbols represent our six regions of integrated 
spectra. 

The V = 1-0 transitions are clearly well thcrmalized in 
all regions, consistent with a similar thermal excitation 
temperature. By a linear least square minimisation to 
the V = 1-0 transitions, we have found a common excita- 
tion temperature of 1930 ± 114K which is in agreement 
with the temperature range 2000 — 3000 K derived from 
Fig-Hni The reason for using only the v = 1-0 lines for 
the fit resides in that these lower energy transitions are 
mainly susceptible to thermal processes, while higher en- 
ergy transitions {v — 2-1) are more easily affected by fluo- 
rescence and other non-thermal mechanisms (if present) . 

For the v — 2-1 transitions, however, the situation 
looks dissimilar. While the (2,5) transitions are consis- 
tent with the thermal linear fit, the (2,3) ones are shifted 
upwards this relation. The exception is region A for 
which all transitions are consistent with a thermal mech- 
anism. The (2,3) deviation of the remaining areas could 
be explained by some contribution from UV-fluorcscence, 
expected to be present in starforming regions. However, 
when fluorescence dominates the H2 excitation, a larger 



The Local LIRG NGC5135 



17 



TABLE 4 

H2 COLUMN DENSITIES DERIVED FOR SIX REGIONS IN NGC 5135. 





Region A 




Region B 




Region C 




Region D 




Region E 




Region F 




Line 


b 


N'^ cor 


b 


COT 


b 


cor 


b 


COT 


b 


cor 


b 


COT 












[1021cm-2j 




[1021cm-2] 












i-U 


Q /1 

o.4y 


y.ou 


111 
1. i i 


i.o / 


U.oo 


l.DU 


1.0/ 


Z.oD 


1 QQ 

i.yo 


0. io 


U.oo 


1.01 




±0.66 


± 2.17 


±0.26 


±0.49 


±0.28 


±0.56 


±0.30 


±0.54 


±0.38 


±0.73 


±0.20 


±0.41 


l-OS(l) 


10.54 


30.38 


2.78 


4.85 


2.45 


4.68 


4.51 


6.97 


6.94 


11.65 


2.37 


4.38 




±0.46 


±4.64 


±0.22 


±0.78 


±0.25 


±0.81 


±0.24 


±1.05 


±0.33 


±1.73 


±0.18 


±0.73 


1-0 S(2) 


3.16 


9.78 


0.90 


1.63 


1.07 


2.13 


1.33 


2.12 


2.41 


4.20 


0.83 


1.60 




±0.40 


± 1.96 


±0.21 


±0.45 


±0.27 


±0.63 


±0.21 


±0.46 


±0.30 


±0.82 


±0.17 


±0.41 


1-0 S(3) 


9.10 


30.20 


2.05 


3.86 


2.11 


4.40 


3.43 


5.62 


6.11 


10.99 


1.93 


3.86 




±0.38 


±5.17 


±0.19 


±0.71 


±0.23 


±0.85 


±0.21 


±0.96 


±0.29 


±1.82 


±0.16 


±0.72 


2-1 S(l) 


0.77 


2.04 


0.40 


0.66 


0.38 


0.68 


0.60 


0.90 


0.92 


1.48 


0.26 


0.46 




±0.25 


± 0.70 


±0.13 


±0.24 


±0.17 


±0.31 


±0.15 


±0.25 


±0.30 


±0.52 


±0.09 


±0.17 


2-1 S(3) 


0.54 


1.64 


0.25 


0.44 


0.15 


0.29 


0.26 


0.42 


0.38 


0.64 


0.16 


0.29 




±0.22 


± 0.70 


±0.11 


±0.21 


±0.12 


±0.24 


±0.10 


±0.17 


±0.15 


±0.27 


±0.08 


±0.16 



(") Total column densities derived by using the observed H2 line fluxes. Below the main values, their ±1 a errors. 

C") Total column densities after extinction correction for each observed H2 line and region as listed in Table[7] Below the main values, their 
±1 a errors. 



positive offset between the 2-1 transitions and our fit 
would be expected, with fine ratios hkc 2-lS(l)/l-0S(l) 
of about 0.5. The observed values in Fig. [TO] are rather 
smaller (< 0.2), suggesting that UV- fluorescence might 
be present but does not dominates the H2 excitation. 
The (2,5) transition case is rather different, with larger 
uncertainties in the data. However we just point out 
that almost all data points lie very close to the thermal 
linear fit, which is not expected for v = 2-1 transitions 
if fiuorescence is present. As mentioned in Sec. 13. 5. 11 
the Draine & Woods (1990) X-ray models predict a (2,5) 
transition decrement which may explain the observed be- 
haviour but, in principle, only for regions A and E (the 
only two with clear X-ray sources). The existence of a 
second thermal component for the v = 2-1 transitions 
cannot be explored given the uncertainties. 

3.6. [Fell] as a SNe tracer 

For many years different authors have suggested that 
the [Fell] -emission from galaxies traces the fast shocks 
produced by SNR and so, their SNe activity (e.g. Graham 
et al. 1987, Moorwood & Oliva 1988, Oliva et al. 1989, 
Greenhouse et al. 1991, Colina 1993, Alonso-Herrero et 
al. 2003, Labrie & Pritchet 2006). The partial ISM ion- 
ization produced by this mechanism is suitable for most 
of the Fe atoms to be in low ionization states. 

Different authors have derived empirical correlations 
we can use to estimate the SNe rate directly from the 
[Fell] data. Calzetti (1997) has found the following ex- 
pression for [Fell] measurements 

SNrate,Ca97{yr-^) = 5-38 • L[p,jj]{W) X lO'^s, (5) 

for which SNes release 10^^ erg of energy during their 
lives. Also, from their study of SNR in M82 and 
NGC 253, Alonso-Herrero et al. (2003) have found that 



SNrate,AHQ3{yr-^) = 8.0 • £[Fe//](W^) X 10" 



(6) 



In Tabled] we present the derived SNe rates for regions 
B-F hy using eqs.jMS] Region A has been excluded from 
the analysis because it is likely that the [Fell] -emission 
comes from the AGN and not from SNe. 



TABLE 5 

SNe Rates derived from [FeII] emission 



Reg 


Ca97 


A-H03 








Bl' 


0.006 ± 0.001 


0.009 ± 0.002 


C* 


0.006 ± 0.001 


0.009 ± 0.002 


D>' 


0.011 ± 0.002 


0.017 ± 0.004 


If 


0.027 ± 0.006 


0.041 ± 0.009 


pi 


0.006 ± 0.001 


0.009 ± 0.002 


Tot. -A'''" 


0.368 ± 0.084 


0.547 ± 0.125 



(") The complete H-band FoV was collapsed in a single spectrum 
and the measured [Fell] flux was extinction corrected by our me- 
dian FoV-extinction value of Ay = 6 mag. For this calculation, 
the [Fell] flux from region A has been subtracted from the total 
because it likely comes from the AGN and not from SNe. ('') Our 
extinction corrected [Fell] fluxes were used. For B-F it corresponds 
to apertures of 157 pc in diameter. 

We remind the reader that all SNe rates derived from 
[Fell] are upper limits because we are assuming that all 
[Fell] flux comes from SNR. We have already seen in 
previous sections that there exist other processes that 
might contribute to [Fell] enhancement. 

The total-FoV SNe rate from TableO is rather high 
compared to other well studied starburst galaxies. The 
rates derived for the entire galaxies M82 and NGC 253 
by Alonso-Hererro et al. (2003) of « 0.1 yr~^ are about 
5 X lower than our [Fell] estimation between the central 

2.3 kpc X 2.3 kpc of NGC 5135. The source of this 
discrepancy is the different SFRs of these galaxies. 
By using Sanders et al. (2003) IRAS data and the 
Kennicutt (1998) law, we estimate total SFRs of 15.0, 

3.4 and 7.O]Vl0yr-i for NGC 5135, NGC 253 and 1VI82, 
respectively. For NGC 5135, we subtract the AGN 
contribution (25%, Alonso-Herrero et al. 2006b) from 
the total far-IR luminosity used in this calculation. As 
we can see, the different SNe rates between NGC 5135 
and NGC 253 can be explained by the stronger SFR of 
the former. The relative SFRs of NGC 5135 and M82 
can also explain their different SNe rates between a 
factor of 2. Such a difference can be easily attributed to 
the many uncertainties of this calculation, including the 
flux measurements, the presence of different alternative 
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Fig. 12. — [Fell] flux map. The contours are the 6 cm VLA 
radio emission from Ulvestad & Wilson (1989). 

sources for those fluxes and the implicit models assump- 
tions for the SNe rates and SFRs determination. 

Forbes & Ward (1993) have found a tight correlation 
between the 6 cm and [FeII]-emission and it is also well 
known that SNe and SNR shells are strong radio emit- 
ters. Therefore, to first order, a spatial correspondence 
between radio and [Fell] -emission could be expected. To 
test this, in Fig. [T^ we compare our [Fell] flux map with 
VLA 6 cm radio emission (contours, 0.91" x 0.60" beam 
resolution) obtained by Ulvestad & Wilson (1989). Both 
images were aligned by matching the corresponding nu- 
clear {A region) unresolved peak. As we can see, the spa- 
tial structures are in good agreement, suggesting that at 
least part of the [FeII]-emission comes from SNRs. 

We can go further on our comparison by contrasting 
the SNe rates derived from radio and [Fell] -emission in 
our B-F regions. To do so, however, we have the limi- 
tation of the poorer spatial resolution of the 6 cm VLA 
image. Therefore, we decided to keep the centers of our 
B-F regions but use a larger circular aperture of 0.91" 
(« 260 pc in diameter, the larger dimension of the radio 
beam size). We label our 0.91" aperture regions as B'-F. 
In this way, both our [Fell] and 6 cm emission were re- 
measured within the new apertures. We highlight that, 
for this larger aperture, regions 77 and share about 
40 % of their spaxels. 

In a similar way to [Fell], other works have explored 
the relation between the 6 cm radio flux and the SNe rate. 
By using models and empirical correlations, the different 
authors have made diverse assumptions concerning the 
origin of their radio observations. From their study of 
compact radio sources in M82, Huang et al. (1994) de- 
rived the following expression 

SN^:Z,Hu9iiyr-') = 0.22 • L^iWHz-') x IQ-^^. (7) 
These authors assume that the entire non-thermal radio 



emission of M82 comes from SNRs. L^, is the 6 cm radio 
flux in [WHz-i]. 

In a study of the brightest radio sources in Arp 299, 
Neff et al. (2004) derived a similar expression for the 
SNe rate as a function of the 6 cm radio flux from SNRs 

According to the authors, this equation is valid for few 
hundred parsec spatial scales, the same order of our own 
aperture diameters. Eq.[8] has been modified to adapt 
the measured radio flux in NGC 5135 to the distance of 
Arp 299 (« 42 Mpc, for which the original expression was 
made). The spectral index a was chosen to be —0.58, a 
value consistent with typical Galactic SNRs (Weiler et al. 
2002) and with the Neff et al. measurements in Arp 299. 
The parameter / is the fraction of radio flux density that 
can be associated with individual SNRs on parsec scales. 
Because there exist other sources of radio emission apart 
of SNR, we adopt a / value of 0.15 following Neff et 
al. An / value of 0.1-0.2 provides a calibration similar to 
those obtained by Condon (1992) and Huang et al. (1994) 
works on radio sources in starburst galaxies. Also, the 
chosen value of / is in excellent agreement with Colina 
& Perez-Olea (1995) models: for the majority of them 
the fraction of SNR radio emission with respect to other 
sources (free-free emission, young radio SNe) ranges be- 
tween 9-19% of the total flux. 

Finally, we have also used Perez-Olea & Colina (1995) 
model predictions for the SNe rate, based on our 6cm 
measurements and Br7 age estimations. An instanta- 
neous burst and Salpeter IMF have been assumed, and 
different masses M„p (30, 60, 1200 Mq) and Mjou; (0.85, 
3Mq) have been considered for the calculations. These 
models make a detailed consideration of the different ra- 
dio components, including free-free radiation, young ra- 
dio SNe and older SNRs. In Table[6]we present the SNe 
rates derived from [Fell] and 6 cm radio emission (eqs. [SI- 
IS]) for regions B'-F . 

Table[6] clearly shows the excellent agreement between 
[Fell] and radio SNe rate predictions for these ~ 200 pc- 
scale regions. This agreement is totally independent of 
the correlation or model used to estimate the SNe rate. 
These results provide additional support to previous find- 
ings in normal starburst galaxies like M82 and NGC 253 
where similar SNe rates are predicted from [Fell] and 
6 cm radio emission. The radio predictions in this table 
are of the same order as those of the brightest radio re- 
gion in Arp 299 (0.5-1.0yr-\ Neff et al. 2004), a bright 
IR luminous galaxy. 

4. SUMMARY 

In this paper we have studied the central 2.3 kpc 
of NGC 5135, a nearby LIRG galaxy. By using new 
NIR IFU spectroscopy (VLT/SINFONI) we have traced 
different structural components including the coronal 
gas (an AGN-activity tracer), ionized (young stellar 
population) and partially ionized gas (SNe remnants), 
warm molecular gas and the giant/supergiant stellar 
population. We have selected individual regions of 
interest for a detailed analysis (see Fig.[3]top panels, and 
Fig. m top- left panel, for their spatial identification). By 
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TABLE 6 

SNb Rates derived from [FeII] and Radio emission 



Reg" 


Ca97[Fen] 


A-H03[FeII] 


Hu946cm 


Nc046cm 






[yr-^] 










B' 


0.012 ± 0.002 


0.017 ± 0.004 


0.012 ± 0.004 


0.019 ± 0.007 


0.013-0.014 


a 


0.012 ± 0.003 


0.018 ± 0.004 


0.010 ± 0.005 


0.016 ± 0.008 


0.011-0.012 


D' 


0.024 ± 0.005 


0.036 ± 0.008 


0.032 ± 0.075 


0.051 ± 0.118 


0.036-0.037 


Ef 


0.051 ± 0.012 


0.076 ± 0.016 


0.064 ± 0.001 


0.100 ± 0.002 


0.071-0.074 


P 


0.011 ± 0.002 


0.016 ± 0.004 


0.012 ± 0.004 


0.019 ± 0.006 


0.018-0.019 



(") Same centers as B-F regions with aperture of 0.91". ('') For Perez-Olea & Colina (1995) models, we present ranges of SNe rates for 

each region, where different upper- and lower-limit masses have been considered for a Salpeter IMF. 



using these data, we have shown the potential of the IFS 
technic when applied to a prototype LIRG: studying the 
complex spatial structure of its different components, the 
different gas ionization/excitation mechanisms, the AGN 

activity, the SNc rates, the internal gas/dust extinc- 
tion and the stellar population ages. Our conclusions are: 



for a galaxy with recent star formation, UV-proccsses 
like fluorescence have, at the most, a rather marginal 
contribution even in Hli regions {B, C and D). The 
column densities found (Nh ~ 10^^-10^^ cm~^) are 
consistent with previous literature estimations at other 
wavelengths. 



1) Overall, we have found that SNR shocks play a 
dominant role on ionizing/exciting the different gas 
phases in the central 2.3 kpc of NGC 5135. Our NIR 
data reveals that other processes associated to recent 
star formation and AGN activity usually have a rather 
localised and smaller overall impact. By using Br7- 
emission, we estimate the AGN contribution to the gas 
ionization in < 3% of the total within our FoV. In the 
same area, the AGN contributes in < 3-5% to the total 
NIR luminosity. 

2) We have found the largest (~ 600 pc in projection) 

[SiVI]-ionization cone reported in literature until now. 
This structure is at least 4x larger than any previous 
detection in active galaxies. Pointing in opposite direc- 
tion, a fainter counter- [SiVI]-cone has been also detected. 

3) The SNe rates derived from the [Fell] emission 
are in excellent agreement with 6 cm radio emission 
predictions, reinforcing the use of [Fell] as a SNe activity 
tracer. Typical rates between 0.01-0.04 yr~^ were found 
for individual ~ 200 pc regions. The [Fell] prediction for 
our entire FoV (0.4-0.5 yr"-'^) is rather high compared to 
other starburst galaxies. Such a difference, however, can 
be explained by the higher SFR of NGC 5135. 

4) In a complex spatial distribution, different pro- 
cesses drive the gas ionization depending on the region 
of interest. These include photoionization from star- 
forming knots (Hii regions 5, C and D), SNc remnant 
shocks {E and F) and AGN photoionization (region A). 
Mechanisms like stellar X-rays photoionization and star 
formation might also be present in other regions {E and 
A, respectively). 

5) The molecular gas excitation mechanisms are 
mainly thermal (~ 1930 ± 114 K) in all galaxy regions, 
being X-rays and SNe remnant shocks the dominant 
processes. Contrary to the expectations we might have 



6) As upper limit estimations, the last burst of star 
formation occurs about 6-8 Myr ago in all regions, 
consistent with previous findings in the UV/optical. 
The presence of a second, older stellar population dom- 
inated by bright red giant /supergiant stars is strongly 
suggested by the data. Although NIR-continuum, 
COA2.29/xm and NaIA2.21/xm observations agree with 
this idea, the simple stellar population models cannot 
clearly discriminate between both populations. 

7) For the nucleus (region A), the presence of deep 
COA2.29/zm and NaIA2.21/zm features highlight the 

importance of star formation in the AGN neighbourhood 
(< 90 pc from the center). This is also consistent with 
previous findings in X-rays. 

8) A mean internal visual extinction of Ay = 6 mag 
(« 0.55 mag at 2.2 /im) has been measured in the central 
2.3 kpc of this galaxy. For the nucleus (region A) a higher 
mean value of Ay = 10 mag was found, with a peak of 
17 mag in the very center of this area. This result is 
qualitatively consistent with high extinctions reported 
for this region in the literature. 
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Fluxes, FWHM and EW for six regions in NGC5135 
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^ Total measured fluxes within aperture of 0.62" (0.54" for [Fell] in ff-band) line fluxes. 

^ Galactic extinction has not been applied given the low value of E(B — V) = 0.06 (Schlegel et al. 1998). Total fluxes corrected by local extinction {{Ay), see Sec l3.2l l. 
" Stellar absorption features. CO band corresponding to the 2-1 transition at 2.29 fj,m. NaT doublet at 2.206-2.208 /^m. 



